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Abstract 

The new experiment planned at Brookhaven to measure the anomalous magnetic mo¬ 
ment of the muon a^ = ( g^ — 2)/2 will improve the present accuracy of 7 ppm by about a 
factor of 20. This requires a careful reconsideration of the theoretical uncertainties of the 
g — 2 predictions, which are dominated by the error of the contribution from the light quarks 
to the photon vacuum polarization. This issue is crucial also for the precise determination 
of the running fine structure constant at the Z-peak as LEP/SLC experiments continue to 
increase their precision. In this paper we present an updated analysis of the hadronic vac¬ 
uum polarization using all presently available e + e _ data. This seems to be justified because 
previous work on the subject was based to some extent on preliminary or incomplete experi¬ 
mental data. Contributions from different energy ranges are presented separately for g — 2 of 
the muon and the r-lepton and for a(M§). We obtain the results ajj ad * = (725±16) x 10 -1 ° 
and ai) ad * = (351 ± 10) x 10 8 , where the asterisk indicates the dressed (renormalization 
group improved) value. For the effective fine structure constant at Mz = 91.1888 GeV we 
obtain Aaff d = 0.0280 ± 0.0007 and a(M§) _1 = 128.896 ± 0.090. Further improvement in 
the accuracy of theoretical predictions which depend on the hadronic vacuum polarization 
requires more precise measurements of e + e _ cross-sections at energies below about 12 GeV 
in future experiments. 


1. Introduction 


The anomalous magnetic moment of the muon has been measured with very high precision at the CERN Muon 
Storage Ring (ij. R is one of the best measured quantities in physics. As it can be calculated with high accuracy 
[^, [| [|] it provides an extremely clean test of electroweak theory and may give us important hints on possible 
deviations from the Standard Model (SM)]|, ||, ||]. The special interest comes about because at {t = e, ft, r) 
corresponds to a helicity flip coupling £^a £r which must vanish at the tree level for any fermion in any 
renormalizable field theory. In the SM it is thus a finite calculable quantity nonvanishing only due to quantum 
fluctuations. Another interesting feature is that at is finite only if the WWq-coupling is of the Yang-Mills type. 
In the limit mt —> 0 regular contributions to at vanish in the SM because £l does not couple to £r in this limit. 
This brings one power of mt into the effective U') vertex; a second power is due to the normalization to the 
lepton Bohr magneton e/2 mt- 

The mf-dependence of the weak interaction and the vacuum polarization effects makes them completely unob¬ 
servable for the electron. The effects are enhanced however by the large factor for the muon relative to 

the electron. This enhancement is very welcome and happens for all kinds of new physics that could couple to 
photons and leptons. Therefore a„ is an important tool for obtaining stringent upper bounds on new physics 
contributions. For the r (see e.g. g], [l(J) there is an additional enhancement factor m%./mi which magnifies the 
interesting effects. However, because of the short r-lifetime, a measurement of a T is very difficult. Therefore 
one is still off by a few orders of magnitude in establishing its value experimentally. The best possibility to 
determine a T is by the leptonic radiative r-decay t ~ —> e~v e v T 7 which also provides the best present upper 
bound fnl . 

The present accuracy Q 

a axp = (11659230 ± 85) x 10” 10 

allows us to test QED with very high precision. The theoretical prediction includes a number of terms of 
different origin. We may write 


,,the _ QED had , weak . new 

where the first, and by far largest term, is from pure QED, a kad denotes the virtual hadronic (quark) contri¬ 
bution, which will be studied in this paper, a™ eak summarizes the SM effects due to virtual W, Z and Higgs 
particle exchanges and a aew stands for possible contributions from extensions of the SM. If we assume the last 
term to be zero we have the theoretical prediction j^] 


a/ he = (11659192 ± 18) x 10“ 10 , 


and therefore at present we have 

a exp _ a the = ( 38 ± 87) x 10 -10 _ 


With the precision attempted by the forthcoming Brookhaven experiment 112 one should be able to establish the 
weak SM contribution a™ eak ~20xl0~ 10 , for example. In fact, however, such contributions may be concealed 
by the theoretical uncertainty which is of a similar size. The latter is dominated by the uncertainty of the 
hadronic contribution and therefore a careful analysis of the problem is of primary importance. The problem is 
that the low energy hadronic effects cannot be calculated in perturbative QCD and one has to rely on the semi- 
phenomenological dispersion theoretical approach |jl3], 14 |n| which allows us to compute a kad as an integral 
over experimental data from e + e _ annihilation, 
theoretical uncertainty of g — 2 predictions. 


The experimental errors of the data of course imply then a 


A comparison of different results based at least partially on this approach is given in Tab. 1 which illustrates 
the present status. 
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Table 1: Comparison of estimates of a dad • 10 10 by different authors 


a ha d . 1Q 10 

Author 

663 ± 85 

Barger et al. 

684 ± 11 

Barkov et al. 

707 ± 6 ± 16 

Kinoshita et al. 

710 ± 10 ± 5 

Casas et al. 

705 ± 6 ± 5 

Martinovic, Dubnicka 

724 ± 7 ± 26 

Jegerlehner 

699 ± 4 ± 2 

Dubnickova et al. 

702 ± 6 ± 14 

This work 

725 ± 6 ± 15 

RG improved 


Year [Reference] 


1975 

H 

1985 

0 

1985 

0 

1985 

@ 

1990 

@ 

1991 

0 

1992 

[0 


The second quantity considered in this paper is the hadronic contribution Aa{^ d (M§) of the 5 light quark flavors 
to the shift of the fine structure constant from the Thomson limit to the Z-resonance Aa = 1 — a/a(M§). 
This quantity can also only be estimated reliably in terms of the experimental e + e _ data with corresponding 
uncertainties. The effective fine structure constant a(s) plays a crucial role in precision physics at all energy 
scales beyond the very low energy region. It is particularly important for all physics of the weak gauge bosons, 
as being studied currently at LEP and SLC. Almost all SM predictions of observables in terms of a , G^ and Mz, 
the most precise set of input parameters available, depend on Aa(s). Again, the uncertainty of Aa|^ d (Mf) 
is a limiting factor which obscures the interpretation of precision measurements at a certain level. In fact, 
for the leptonic effective weak mixing parameter sin 2 0^ t (Mz), which is determined in e + e _ —> l + l~ at the 
Z-resonance, the LEP collaboration |23j has reached an experimental error <5 sin 2 l/~ 0.0004 which is only 
slightly larger than the hadronic uncertainty of about 0.0003 of the SM prediction for this quantity. It thus has 
become a very crucial question whether it is possible to reduce the error of Aa|^ d (Mf) further and what the 
perspectives are in the future. The present situation is summarized by the results in Tab. 2. 

Note that some of the values have been shifted by 

|-f (1 + ^) In (M' z /Mz) 

from the pre-LEP reference value Mz = 93 GeV to the current value of Mz- Ref. |2f| quotes Aa[ 5 ^(— Qq) = 
0.0145 ± 0.0012 for Qq = 79 GeV 2 . We have added Aa[^ d (M§) — Aaj^f— Qq) = 0.0138 which is obtained by 
using 3rd order perturbative QCD. For a comparison of the earlier results |lo) we refer to [Q. The differences 
are mainly due to a different treatment of the systematic errors and/or different model assumptions. The present 
situation justifies a reconsideration of the problem. 


Table 2: Comparison of estimates of Aa^\(M'^) by different authors 


A4ad(^D 

Author Year [Reference] 

0.0285 ± 0.0007 

0.0283 ± 0.0012 

0.0287 ± 0.0009 

0.0282 ± 0.0009 

0.02666 ± 0.00075 

0.02732 ± 0.00042 

0.0280 ± 0.0007 

Jegerlehner 1986 (tJ 

Lynn et al. 1987 j2f| 

Burkhardt et al. 1989 |2Cj] 

Jegerlehner 1991 j2?j] 

Swartz 1994 |2§|] 

Martin and Zeppenfeld 1994 ]2S|] 

This work 
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To calculate the necessary dispersion integrals below we prefer to use direct integration over the experimental 
values of cross sections. In this approach one can take into account uncertainties of separate measurements 
in a straightforward manner. The alternative method which was used in most of the previous works is to 
make a fit of the experimental points within some model and integrate the arising parametrization of the data. 
This procedure inevitably leads to a model dependence and it is not clear how experimental errors especially 
systematic uncertainties can be taken into account. 

The two recent papers |2^, ^9| appeared while we were writing up this update. We have added some remarks 
which should clarify at least part of the discrepancies between these papers and the present analysis which 
yields results consistent with previous ones (e.g. |2?|] ). 

In Secs. 2 and 3 we will briefly review how the hadronic contributions to an and Aa are determined by the 
e + e _ data. Details on the evaluation of the dispersion integrals and the treatment of the errors are discussed 
in Sec. 4. In Sec. 5 we describe the e + e _ data which we use in our analysis and in Sec. 6 we present the results. 
A brief summary and outlook follows in Sec. 7. 

2. Hadronic contributions to the anomalous magnetic 
moments of the leptons 

The leading hadronic contribution to g — 2 is due to the photon vacuum polarization insertion into the vertex 
diagram of the electromagnetic vertex of a lepton. The corresponding diagram is shown in Fig. 1. 



Fig. 1: Leading hadronic vacuum polarization contribution to g — 2 of a 

lepton. 


The “blob” represents the irreducible photon self-energy. Subleading hadronic contributions are obtained by 
multiple self-energy insertions or insertions of irreducible light-by-light scattering (4y) or higher amplitudes 
ii- We will not discuss such subleading terms in this paper unless stated otherwise. 

As already mentioned before, the low energy hadron effects which are needed here, and in principle are de¬ 
termined by QCD, cannot be obtained by using perturbation theory. Fortunately, this contribution may be 
calculated in terms of the experimental total cross section cthad = cr(e + e _ —> hadrons) of e + e _ annihilation into 
any hadronic state by using the familiar dispersion integral |L3|, [m|, |l^| 


a 


had 


l 

47T 3 



ds a hld( s ) K ( s ) 



R(s) K(s) 


(1) 


which can be evaluated by using the experimental data for (Jhad(s ) or R(s) up to some sufficiently high energy, 
e.g., E cut = 40 GeV, and by perturbative QCD for the high energy tail. 
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For the moment we are interested in calculating the contribution from the irreducible photon self-energy, in 
which case we have to use the “undressed” (= lowest order with respect to QED) hadronic cross section 


<r { had( s ) = Vhadis) (a/a(s)) 2 


( 2 ) 


We refer to the Appendix for a brief discussion of “dressed” versus “undressed” quantities in dispersion relations. 

The second form in Eq. (;lj) is convenient for the evaluation of contributions at higher energies. It is an expression 
in terms of the cross-section ratio 


R(s) = 


0 ’tot(e + e —> 7 * —> hadrons ) 




■+ P - 


H+H ) 


(3) 


Note that R(s), by the proper definition according to Eq. is the ratio of the total cross sections which is 
determined by QCD. In perturbative QCD we have 132 


R(s) = 3 Q 2 f\j 1 ~ 4mf /s (1 + 2 m 2 /s) (l + a + c\a 2 + C 2 a 3 + • • •) 


(4) 


where Qf and to/ denote the charge and mass of the quark, respectively, a = a s (s)/ir with a s (s) the strong 
coupling constant, and 


ci = 1.9857 — 0.1153AI/ 

c 2 = —6.6368 — 1.20027V/ — 0.0052JV 2 — 1.2395 (E«/) 2 /(3£«?) 


in the MS scheme. Nf is the number of active flavors. This result is applicable about 1 GeV above the 
resonances and at sufficiently high energies and will be used in particular to calculate the high energy tail of 
Eq. (jlj). We will assume a top mass of 173 GeV Q for the calculation of the perturbative tail. As a check 
we also will compare the experimental data in the region above the ^-resonances with the calculation of R(s) 
presented recently in Ref. 133 , which is improved by charm and bottom mass effects. 

Usually experiments do not determine R as a ratio of the total cross sections as given by Eq. (pj). Rather 
the hadronic experimental cross section is first corrected for QED effects |f4[ [j(| [57J, which include 

bremsstrahlung as well as vacuum polarization corrections. The latter account for the running of the fine 
structure constant a(s). After these corrections have been applied a tot. is divided by the Born cross section 


0o (e 


+ /3 _ 


— \ 47ra: 

3s 


h + h ) = 


so that 


R(s) = 


CTtot(t 


7 * —> hadrons) 


corr 

exp 


ao(e + e —> 7 * 


H+H ) 


Note that, the experimental cross section cr(e + e _ —^> 7 * —> n + ^~) never appears here and is used by careful 
groups to check how good normalization is (see e.g., JH]). The question of how R has been determined precisely 
in a given experiment is not always clear. We will comment on this point when discussing the data below. 

A renormalization group (RG) improvement of the result may be obtained by resumming the multiple irreducible 
self-energy insertions, which is equivalent to using the full photon propagator in Fig. 1. As shown in the 
Appendix, this is simply achieved by using the physical cross section <Jhad(s) or, equivalently, the “dressed” 
i?-function i?(s) dressed = R(s)(a(s)/a) 2 under the dispersion integral Eq. (lj). 

Turning back to Eq. ®): the kernel I\ (s) may conveniently be written in terms of the variable 


_ l-fl, 
1 + P/J. 


Pn = yjl-4rnl/s 


and is given by 


A^) = y (2-:r 2 ) + 


(1 + £ 2 )(1 + x ) 2 


ln(l + x) - x + — ) + 


(4 + x ) 

{1-x) 


x 2 In (a;) 


(5) 
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Note that the function 


- 3s 

K(s) = —^K(s) 


is bounded: it increases monotonically from 0.63 at threshold s = 4m z to 1 at s —> oo. It should be noted that 
for small x the calculation of the function /\ (s), in the form given above, is numerically instable and we instead 
use the asymptotic expansion (used typically for x < 0.0006) 


K{s) = { - + 


17 

12 


+ 


11 

30 


+ I + -^x | x | x | x | x + 


1 + X 
\ — X 


x 2 ln(a;) . 


Other representations of K(s), like the simpler-looking form 


1 


K(s) =-r H —r (r — 2) ln(r) + I 1 — 2r H— r z ) In(x)//? y 


1 


with r = s/mji, are much less suitable for numerical evaluation because of much more severe numerical cancel¬ 
lation (even less stable is the representation utilized in f22f). 

The representation Eq. (j^) of /\(s) is valid for the muon (or electron) where we have s > 4ni 2 in the domain 
of integration s > 4mJ, and x is real, and 0 < x < 1. For the r Eq. (|j) applies for s > 4 m~.. In the region 
4m^ < s < 4where 0 < r = s/m z < 4, we may use the form 


K(s) = — — r + — r (r — 2) ln(r) — ( 1 — 2r + 


tp/i 


( 6 ) 


with w = \J4jr — 1 and <p = 2 tan 1 (w). 


3. Hadronic contributions to the running of the effective 
fine structure constant 


The effective fine structure constant at scale is given by 


a(s) 


a 

1 — A a(s) 


(7) 


where a is the fine structure constant and Aa is the photon vacuum polarization contribution. In terms of the 
one particle irreducible photon self-energy ILy(s) = sll' (s) we have 

Aa(s) =n;(0)-Ren;(s) , (8) 

which for s = M |, for example, is large due to the large change in scale going from zero momentum (Thomson 
limit) to the Z-mass scale /.i = Mz- In perturbation theory, the leading light fermion (m/ -C Mw,\fs) 
contribution is given by 

Aa(s) = 

f 
a 
3n 
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with Qf the fermion charge and N c f the color factor, 1 for leptons and 3 for quarks. We distinguish the 
contributions from the leptons, for which Eq. (|9j) is appropriate, the five light quarks and the top 

Aa = A ai + Aa^ + Aa top ■ (10) 

Since the top quark is heavy we cannot use the light fermion approximation for it. A very heavy top in fact 
decouples like 


A (Xt op — 


a 4 M§ 
3-7T15 m 2 


0 


when mt 3> Mz- 

A serious problem is the low energy contributions of the five light quarks u, d, s, c and b which cannot be 

(5) 

reliably calculated by using perturbative QCD. Fortunately, again one can evaluate this hadronic term Aa^ 
from hadronic e“*"e~-annihilation data by using a dispersion relation together with the optical theorem which 
results in the integral |3f], |3(J 


A^ad (M 2 Z ) = 


Ml 

47 r 2 a 


Re 


ds 


7 (0) 

had 


W 


/4 ml 


s — Ml — ie 


-tr Re 


R(s) 


Kml dS s(s-M 2 Z -ie) 


( 11 ) 


which is very similar to the one we encountered for g — 2 in Eq. (jl|). The only difference is the different weight- 
function multiplying R(s) under the integral. Since Eq. (jlj) has an extra factor 1/s at low s the low energy data 
play a dominant role for while Eq. (11) gets significant contributions from a broad energy range up to about 
E ~ Mz/ 2, as we shall see below. 

Note that the remarks made earlier about the proper definition of R(s) apply here as well. 

Above the T energy we apply a correction factor for 7 — Z mixing to R(s) exp . The correction for the Z-exchange 
contribution is given by the ratio cqq j c^z where 


C QQ 


= ^ <2/ = 11/9 ; c 7 z = cqq - 2 v e CQ v P(s) + c e cfP(s ) 2 


/ 


We denoted the Z -pole factor by 


p (r ) — y/2G M Ml s 

167TQ S — M 2 Z 

and the coefficients are determined by the Zf f couplings Vf and a/ as follows: 


CQv 

= 

E / Qf v f 

= 7/3 -44/9 sin 2 <d w 

Ce 

= 

vl + a 2 e 

= 2 — 8 sin 2 Qw + 16 sin 4 ©w 

C f 

= 

E M + a 2 ) 

= 10 -56/3 sin 2 Q w + 176/9 sin 4 Q w . 

In our normalization v ( 
sin 2 9 l e ept {M z ) = 0.2322 

11 

— 1 + 4sin 2 0w- We use the LEP values Mz = 91.1888 GeV and sin 2 ©w 
for numerical estimates. 


4. Evaluation of the dispersion integral 


In order to obtain a conservative estimate of the relevant integrals we try to rely on the experimental data as 
much as possible and integrate directly the data points by joining them by straight lines (trapezoidal rule). In 
the low energy region it is customary to present the 7 r + 7 r _ and p —> 7 r + 7 r _ data in terms of the absolute square 
of the pion form-factor |F,r(s )| 2 which is related to the cross section by 


er(e + e 


■) = 


7 T a 2 f3f 


\FAs)f 


or 


Rn 


B 3 

= ^f\F,{s )\ 2 


( 12 ) 
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where (3^ = (1 — 4m^/s) 1 / 2 is the pion velocity. In contrast to previous works which were using Gounaris-Sakurai 
[ pt)| kind of parametrizations we will rely on the direct integration of the data in order not to obscure the error 
estimates. The definition of the pion form factor is similar in spirit to the one for R(s). In principle, it is defined 
as a ratio of hadronic to leptonic cross sections such that it corresponds to a purely hadronic matrix element. 
For most of the low energy experiments (in particular the Novosibirsk VEPP-2M ones) R was calculated as 
follows |l7]|: first the hadronic cross section was obtained. Radiative corrections were applied which included 
lepton vacuum polarization only. And then the corrected hadronic cross section was divided by the muon cross 
section calculated using the constant low energy a. The more precise accounting of using the running a with the 
inclusion of hadronic vacuum polarization leads to a numerically very small (a^) or negligible (Aa) difference. 
For consistency we apply the factor (1 -f 2Acq(s))(a/a(s)) 2 to the pion form factor to account for this missing 
correction. The remarks made here for 7r + 7r _ apply as well for the K + K~ and KgK l form-factors. 


4.1. Resonances 


A few exceptions from the direct integration are the narrow resonances ui, <j>, the J/%[) family (6 states) and the 
T family (6 states). Here we can safely use the parametrization as Breit-Wigner resonances 


. , 127 T r ee MlTRTis) 

vswi*) - M|f^(s- M^) 2 + M 2 R T 2 {s) 


(13) 


or as a zero width resonance 

19t r 2 

°nw{s) = r ee S(s - M R ) (14) 

using resonance parameters from the Review of Particle Properties fll|| . The s-dependent width is defined by 
the imaginary part of the irreducible self-energy of the resonance propagator as Imn(s) = M R T(s). Sufficiently 
far above the thresholds it behaves like 


r(s) ~ s/M 2 R T r where T R =T(M%) 


(15) 


The procedure of calculating widths for the narrow resonances of the J/if) and T families is described in the 
Review of Particle Properties, e.g., in the 1994 edition, on page 1661. The particle data group (PDG) lists 
“dressed” (i.e. physical) widths (T, F ee , • • •) rather than lowest order ones (r^ 0 ), rie\ ■ ■ •). This convention exists 
since the 1988 edition of the Review of Particle Properties |l2), and was triggered by the articles Refs. |t|, 
[ p4| and [p5f which pointed out that, in the past, different experiments have been using different conventions 
which depended on the different treatment of the radiative corrections. A discussion of the radiative corrections 
which were performed by different experiments for the T resonances may be found in Ref. [p3| where consistent 
world averages for the dressed widths are calculated by appropriate rescaling of the peaks of the resonances. In 
Ref. [p5| the J/0 and the T resonances were refitted using state of the art calculations for resonance line-shapes 
as they have been developed for precision physics at the Z-resonance, using Eq. ([L3|) with r(s) = s/M R T/}. 
Note that the convention mentioned above is employed by the authors themselves. The PDG only lists the 
results of original papers and averages them. 

In terms of the physical (“dressed”) widths resonance contributions to the R-values are given by 


Rbw 


9 s P ee M'^T r T{s) 
a 2 (s) M 2 R Ta (s — Ml) 2 + M 2 r 2 (s) 


(16) 


and correspondingly for the narrow resonance approximation. The latter is used only for calculating the reso¬ 
nance contributions to an- It should be noted that the Breit-Wigner formula itself, valid in the vicinity of the 
resonance, is a result of summing over all quark vacuum polarization diagrams with any number of loops. 

For the to and the (j> we proceed as described in Ref. Q (see also Q) and use the relativistic Breit-Wigner 
form with a s-dependent width 


r w (a) 


r(u> —> 3 - 7 T, s) + r(w —> 7 r° 7 , s) + r(w —> 271-, s) 



Br(u) 


37t) 


F3it(s) 

f 3 ,(m 2 ) 
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r*(a) 


+Br( “ A) l§t5§) + Br(u 2,r) l£tsk 

= r((j) —> k + k~ ,s) + r(0 —> k s k l , s) + r(0 —> 371 -, s) + r(</> —> 7r°7, s) + r(</ —> 777 , s) 


= j^r*| 


F 


-ffS-ffL 


(«) 


Fk s k l {MI) 


+Br(4> —>■ 37 t) 


J37t(s) 
^3* (M£) 


+Br((j) —> 7 r° 7 ) 


^7T7 (^) 

_F , 7r 7 (M|) 


Br(<f> —> 777 ) 


F 


7)7 


(fL\ 


Frry{M^) j 


( 17 ) 


where Br(V —><■ X) denotes the branching fraction for the channel X and Fy(s) is the phase space function for 
the corresponding channel normalized such that Fx(s) —> const for s —> oo. For the two-body decays h 7 —> P 1 P 2 
we have Fp 1 p 2 (s ) = (1— (7711 + 77 i 2 ) 2 /s) 3 / 2 . The channel V —> 37 r is dominated by F —> p 7 r —> 37 t and this fact is 
used when calculating ^ 3,7 (s) [Q . Before extracting the width, radiative corrections according to Bonneau and 
Martin |^j or Kuraev and Fadin p7| have been performed which include subtracting the electron contribution 
to the vacuum polarization. In these cases the correction to be applied is (1 + 2A a e (s))(a/a(s)) 2 and not the 
full one. 


4.2. The 7r + 7r threshold region 

Experimental data are poor below about 400 MeV because the cross section is suppressed near the threshold. 
Because of the 1/s 2 weight factor for small s in Eq. (jj]) we have to worry whether there could not be a relevant 
contribution missing. Here results from chiral expansion of the pion form factor can be used (see also 

Ref. HU). 

To a good approximation the relevant vector form factor is given by 

F$ hpt ~ 1 + i < r 2 >1 -s + 4 • s 2 (18) 

with < r 2 >y=0.427±0.010 fm 2 and Cy=4.1^tn 2 GeV -4 . The pion charge radius used here was determined 
from the precise spacelike data in Ref. [|49|. We have used the value obtained from fits with a free normalization. 
The error includes the 0.9% systematic error of the data. The crucial point here is that the threshold behavior 
is severely constrained by the chiral structure of QCD via the rather precise data for the pion form factor in the 
spacelike region. The convergence of the momentum expansion can be improved by using the Pade approximants 
rather than the asymptotic expansion itself which is denoted by [0,2] in Pade terminology. 

Fig. 2: The e + e —> 7 r + 7 r data near threshold compared with the pre¬ 

diction of the pion form factor li^s )) 2 for timelike s = E 2 from 
chiral expansion to two-loop order Q. [1,1] and [2,0] denote Pade 
improvements. 

With the information we have on the expansion coefficients we may obtain the [2,0] form 

F$ HPT ~ 1/(1 - ci • s - (c 2 - c 2 ) • s 2 ) , (19) 

with Ci = jjr < r 2 >y and c 2 = Cy\ or the [ 1 , 1 ] form 

F C HpT ~ (1 + (d - c 2 /ci) • s)/(l - (c 2 /d) • s) , ( 20 ) 

which agree with ( |l 8 | ) when expanded in s up to terms of unknown higher orders. The difference obtained from 
the different representations is the “model” error, uncertainties due to missing higher order terms. 

The results are shown in Fig. 2 and provide a good description of the data in the timelike region. 
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4.3. Estimate of the error 


While statistical errors are added in quadrature throughout in our analysis the systematic errors of an experiment 
have to be added linearly. Usually the experiments give systematic errors as a relative systematic uncertainty 
and the systematic error to be added linearly is given by the central value times the relative uncertainty. For 
data from different experiments the combination of the systematic errors is more problematic. If one would add 
systematic errors linearly everywhere, the error would be obviously overestimated since one would not take into 
account the fact that independent experiments have been performed. Since we are interested in the integral over 
the data only, a natural procedure seems to be the following: for a given energy range (scan region) we integrate 
the data points for each individual experiment and then take a weighted mean, based on the quadratically 
combined statistical and systematic error, of the experiments which have been performed in this energy range. 
By doing so we have assumed that different experiments have independent systematic errors, which of course 
often is only partially true 1 !. The problem with this method is that there exist regions where data are sparse 
yet the cross section varies rapidly, like in the p-resonance region. The applicability of the trapezoidal rule is 
then not reliable, but taking other models for the extrapolation introduces another source of systematic errors. 
It was noticed some time ago in Ref. |5lll that fitting data to some function by minimizing y 2 may lead to 
misleading results^ and we insist on avoiding this kind of problems. 

In order to start from a better defined integrand we do better to combine all available data points into a single 
dataset. If we would take just the collection of points as if they were from one experiment we not only would 
get a too pessimistic error estimate but a serious problem could be that scarcely distributed precise data points 
do not get the appropriate weight relative to densely spaced data point with larger errors. What seems to be 
more adequate is to take for each point of the combined set the weighted average of the given point and the 
linearly interpolated points of the other experiments: 

r=— wiRi 


with total error S tot = 1 /y/w , where w = J2i w i an d Wi = 1 /Sf tot . By Si tot = y Sf sta + <5 2 sys we denote the 
combined error of the individual measurements. In addition, to each point a statistical and a systematic error 
is assigned by taking weighted averages of the squared errors: 



There is of course an ambiguity in separating the well-defined combined error into a statistical and a systematic 
one. We may also calculate separately the total error and the statistical one and obtain a systematic error 
<W = VCT — d 2 tQ . Both procedures give very similar results. We also calculate y 2 = w -i (Ri ~ R) 2 and 
compare it with IV — 1, where N is the number of experiments. Whenever S = \/x 2 /(fV — 1) > 1, we scale the 
errors by the factor S, unless there are plausible arguments which allow one to discard inconsistent data points. 


5. The e + e data and the origin of systematic errors 

Some general comments concerning the R determination are in order. In the ideal case one directly identifies 
each possible annihilation channel and measures its cross section. After that R is obtained as a sum of all 
separate contributions. We will call such an approach an “exclusive” one. One should be careful, however, 
while estimating the resulting systematic uncertainty since systematic uncertainties of separate channels may 

1 If there are known common errors, like the normalization errors for experiments performed at the same 
facility, one has to add the common error after averaging. In some cases we correct for possible common errors 
by scaling up the systematic error appropriately. 

2 The problem addressed in Ref. |5^| is that “The best fits to the data which are affected by systematic 
uncertainties on the normalization factor have the tendency to produce curves lower than expected, if the 
covariance matrix of the data points is used in the definition of \ 2 ”. 


9 






contain common parts like, e.g., normalization uncertainties. There may also be some model uncertainties arising 
from the fact that for each specific channel one has to assume some specific mechanism of the production of 
final particles, for example, for the production of four pions, un r, ai7r, pin: etc., while in reality particles 
can be produced by different mechanisms. In a more realistic case identification of separate channels is not 
possible. One observes a certain number of multi-hadronic events of different types, for example, 2 charged 
particles, 2 charged particles plus n photons, 3 charged particles etc., assumes some mechanism of multi-hadron 
production and calculates within it a detection efficiency of observing any configuration of charged particles 
and photons. Minimizing \ 2 one obtains then the cross sections of separate channels and their sum gives a 
total cross section. In old experiments, pion production with an invariant phase space distribution was used as 
a mechanism of particle production (see, e.g. [[78]). Lately the LUND model |52| has been used for detection 
efficiency determination. Here again additional model uncertainties arise. The value of R determined by such 
an “inclusive” method usually yields the total cross section not including production of two-body final states. 
One should also take into account that there may be production channels “invisible” to some experiments, like, 
for example, those having only neutral particles in the final state. Some of them may be accounted for by using 
isospin symmetry. For example, the cross section of the reaction e + e _ —> 7r + 7r _ 37r° is equal to one half of 
that for e + e“ —> 27r + 27r”7r°, independently of the production mechanism. Such corrections although small will 
influence the central value of the calculated integral and will also contribute to the systematic uncertainty. 

Since 1985 when the paper of Kinoshita et al. [^8| was published, a lot of new experimental data on the R 
measurement in e + e _ annihilation into hadrons has been accumulated. Progress in the low energy range was 
mostly due to experiments at VEPP-2M at Novosibirsk where three groups (OLYA, ND and CMD) provided 
independently the information on different exclusive channels of e + e _ annihilation and at DCI at Orsay coming 
from the DM2 experiment. In the Novosibirsk experiments the center of mass (c.m.) energy range from the 
threshold of hadron production up to 1.4 GeV was studied, whereas Orsay experiments covered the c.m. energy 
range from 1.35 up to 2.3 GeV. The high integrated luminosity collected in these experiments allowed them to 
improve considerably their statistical precision. Larger solid angles of the detectors and the use of electromag¬ 
netic calorimeters providing detection of photons with good energy resolution facilitated the identification of 
the large number of exclusive annihilation channels. 


In the Novosibirsk energy range the cross section is dominated by p , to and (f> resonances, resulting in its strong 
energy dependence. Typical multiplicities in multi-pion production do not exceed five, making the total number 
of accessible channels rather small. Besides that, the cross section of the two-body channels is dominating the 
total cross section. All these circumstances make the usual inclusive procedure of R determination from the 
total number of multi-hadronic events of a different type almost meaningless and subject to large uncertainties. 
Therefore experimental efforts were aimed at the selection of exclusive channels followed by the “exclusive” R 
determination by simply adding contributions of different reactions. The energy range from 1.4 up to 2.3 GeV 
systematically studied by Orsay and Frascati groups is much worse understood. It is clear that higher vector 
mesons play an important role, but their precise parameters are not yet known, thus making a phenomenological 
model approach rather ambiguous (see the recent paper [|53|). The cross section of multi-hadronic channels is 
considerably larger than that of two-body channels. However, since the multiplicity and therefore the number of 
possible annihilation channels is still not too big, the approach can be twofold: studies of exclusive channels can 
be complemented by the independent “inclusive” R determination. At larger energies (above approximately 
2 GeV) data on the cross sections of separate channels are missing, only few exclusive channels have been 
measured (usually in the near resonance region) and there are “inclusive” R determinations only. 


The reaction e + e _ —> 7r + 7r“ was studied with high precision by the two groups OLYA and CMD which found 
good agreement between each other and the results of the joint analysis were published in Ref. (T?)]. In the 
CMD experiment 24 points from 360 to 820 MeV have been studied with a systematic uncertainty less than 
2%. OLYA performed scanning of the energy region from 640 to 1400 MeV with a small energy step and had a 
systematic uncertainty from about 4% at the p-meson peak up to 15% at 1400 MeV. Also used were the older 
data near threshold from OLYA |Q, TOF |5{|, NA7 (56) and VEPP-2M 57 as well as the measurements from 
483 to 1096 MeV by DM1 (5f|. The data are shown in Fig. 3. At higher energies this reaction was systematically 
studied by DM2 |)59f . Older results from pn fjo| , BCF fill , and MEA fj|| at Frascati are also included in the 
analysis. 


The reaction e + e —> 7r + 7r tt° was studied by different Orsay and Novosibirsk groups at the u> and </> meson 
(for these data we used the values of the leptonic widths from the Review of Particle Properties (h) to calcu- 
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late analytically the resonance contributions) as well as in the off-resonance region by CMD |G^], ND |Q at 
Novosibirsk and M2N J65| ], M3N Q, DM1 [[38], and DM2 |i9] at Orsay. 


There are two channels in four-pion production: 7r + 7r _ 7r°7r° and 27T+27T - . OLYA jFij, [nj and ND |64|] who 
scanned the energy region from about 640 to 1400 MeV provided results on both, while CMD |72j measured 
the cross section of the latter in 9 points from 1019 up to 1403 MeV. The values of the cross section determined 
by ND are usually higher than those of OLYA in both reaction channels; however, they are within systematic 
uncertainties which are estimated by the authors to be 15% and 10% respectively for ND and 20% in both cases 
for OLYA. CMD claimed a 10% systematic uncertainty and within it agreed with both ND and OLYA. At higher 
energies we used the data from M3N { 


MEA |67|, DM1 |73fl and DM2 


|75|[ . Since we are interested in 
the values of total cross sections only, the mechanism of particle production is in general not important for our 
analysis. However, one should note that the reaction e + e _ —► ujtt 0 plays an important role in the production 
of the 7T + 7r - 27r° final state. Since to has a 8.5% branching ratio for the 7r°7 decay 41 it should be taken into 
account separately since it will be missed in the inclusive analysis. Accordingly we add a un oB[iij — * 7r°7) to the 
total cross section, a contribution which was previously ignored. 


The cross section of the reaction e + e _ —> 27r + 27r _ 7r° was measured at low energy by CMD |72] and at higher 
energy by different groups. The most precise of them are those by DM1 |7(j] and DM2 ]69|| . Its isospin 
partner e + e _ —> 7r + 7r _ 37r° is much worse studied, but as mentioned above, the rigorous isotopic symmetry 
relation requires that its cross section be two times smaller. Since w 7T7 t is a dominant mechanism of five pion 
production, one should make a correction for it similar to that for the reaction e“*"e - —+ ton 0 . 


Six pion production is much less studied. There are three possible final states of which only 37r + 37r _ and 
27r + 27r _ 27r° have been observed. The isotopic symmetry does not give unfortunately the rigorous relations 
between different isotopic partners [ [ttI ; however, one can expect that the corrections for the missing parts are 
small. The existing data are rather controversial. The 772 |78| measurements are typically characterized by 
much higher values of the cross section hardly compatible with the measurements of DM1 JtiJ and DM2 |7^[ . 
Measurements at low energy are also available from CMD j|72| . 

The reaction e + e _ —> was studied by Q and |8l| . 

Production of kaon pairs was studied by OLYA (8^, |8^[ , CMD Q, DM1 Q and DM2 JstJ. 


DM1 1881 and DM2 |74j measured also cross sections of the reactions e + e —> KKn and KKtttt. However, 
all possible combinations have been studied. The corrections for the missing modes should be applied. 


not 


Information on the production of baryons is scarce. Measurements by DM1 |8S|] , DM2 [Q and FENICE pi] 
showed that the cross section is small. However, it should be taken into account if R is determined inclusively. 


Our compilation of f?, obtained following the procedure described in Sec. [L3], is shown in Figs. 4 to 8, together 
with original compilations by the experimental groups themselves. It should be mentioned that part of the 
large systematic errors given by experiments in particular below the J/ip are supposed to account for the lack 
of understanding or performing the radiative corrections. Most experiments have applied radiative corrections 
including at least the vacuum polarization contribution from the electron, which is the dominant contribution at 
low energies [fl~i| . [l6, 37|. The hadronic contribution to the vacuum polarization has been known since the work 
of Berends and Komen [Q in 1976 but usually it was not included in the QED corrections. We shall assume 
that on the average experiments have subtracted only the electron contribution to the vacuum polarization and 
accordingly rescale the R-values by (1 + 2A a e ){a/a(s)) 2 below the J/ip. As we mentioned before this remark 
applies to the uj and <p resonances as well. 


In the high energy region we distinguish the J/ip and the T resonances and the background inclusive mea¬ 
surements of the total hadronic cross section which is usually presented in terms of i?-values. The resonance 
contributions are taken into account as explained in Sec. 4.1. Masses, widths and the electronic branching frac¬ 
tions are taken from the Review of Particle Properties jll|. Since only total errors are given for these quantities, 
we treat the error from the mass (negligible) and width quadratically and the one from the branching fraction 
linearly within a family. In this way we take into account that the systematic errors should be added linearly. 

In the region from the J/ip to the T ^-measurements are available from Mark I |)5], |96|, DASP |)7]], PLUTO ]98| ], 
LENA j99j, Crystal Ball (CB) |10C[ | and MD-1 [101] (see Fig. 7). The Crystal Ball Collaboration has carefully 
reanalyzed their old data and obtained R(s) values substantially lower than Mark I and in agreement with 
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the other experiments PLUTO, LENA and MD-1. The results are now much closer to expectations from 
perturbative QCD (i? ~ 10/3 at lowest order). The change of the data is mainly due to a up-to-date treatment 
of the QED radiative corrections and r subtraction. Ref. 100 gives a detailed description of how precisely the 
radiative corrections have been performed. The procedure is based on the calculation of Ref. j3(J which has 
been applied by all DESY experiments. Thus, except from the Mark I data, all data above the J/ip resonance 
have published properly normalized R-values. The Mark I values we rescale by (1 + 2A ai)(a/a(s)) 2 , because 
we assume that hadronic vacuum polarization effects have not been subtracted, while the leptonic correction 
has been taken into account in a non-resummed form. 

Below we will study the option of including and discarding the Mark I data in the overlapping range, where the 
Mark I data systematically lie 28% higher as we can see in Fig. 7. 

Above the T we include PETRA and PEP data up to E cut = 40 GeV. In this range we correct for the 7 — Z 
mixing contribution as described earlier. For larger energies the 7 — Z mixing would be substantial and make 
the analysis less transparent. In fact perturbative QCD is reliable for evaluating the high energy contribution 
already at lower energies. In the PETRA/PEP range there are many measurements of rather high accuracy 
available, a fact which is particularly important for the precise determination of a(Mz)- Data are mainly from 
CELLO |102|, JADE |l 03| , MARK J [ 1041 and TASSO |105[. Further included are the ^-measurements from 
DASP f97|7DHHM pcfTCLEO |l07| , CUSB p 8 | , HRS |0§, MAC JTTt 


and MD-1 [1011(see Fig. 8 ). 


6. Results 


The results presented here have been obtained following the procedure described in Sec. 4.3. For the 7 r + 7 r - data 
we have three sets of data points as collected in Ref. |l7j. They are displayed in Fig. 3 and we label them by 
DM1 + (which includes data from DM1, NA7 and TOF), CMD and OLYA. In the range 0.81-1.4 GeV we use a 
compilation similar to the one in Ref. p4j but with all available data which are shown in Fig. 4. This is don e by 


adding up the individual channels, where for each channel weighted averages are taken as described in Sec.4.3 


The (/-resonance is taken into account in analytic form in the narrow interval between 1.00 and 1.04 GeV. 
Outside this interval the <p contribution is included in the background. The cut at 1.4 GeV is justified as it is 
the energy limit of VEPP-2M. From 1.4 to 2.3 GeV we combine R-values published by the experiments MEA, 
772 , BB, M3N, DM1 and DM2; from 2.3 to 3.1 GeV data collected by 772 , DM2 and Mark I are available, where 
the last two give R while the others R(n > 2), in which case the 2-body channels have to be taken into account 
separately. In taking weighted averages here we did not account yet for the fact that some measurements have 
been performed at the same machine and hence have common normalization errors. We therefore enlarge the 
systematic error in this domain by a factor \/2 to be on the conservative side. The data for this region are 
shown in Figs. 5 and 6 . In the region between the J/ip and T resonances we again compare two methods. First 
we combine data by calculating the weighted averages as described above. The Mark I data are treated as an 
independent set, and are eventually combined after integration. Alternatively, we split the region in such a way 
that results from different experiments can be combined after integration. From 3.1 to 3.6 GeV we have data 
from Mark I, form 3.6 to 5.2 GeV Mark I, DASP and PLUTO, from 5.2 to 7.2 GeV Mark I, CB and PLUTO 
and from 7.2 to 9.46 GeV PLUTO and MD-1. Here we take weighted averages after integration. The different 
energy ranges mentioned are treated as independent when adding up the results, a procedure which again may 
not be fully justified. 


The data above the J/ip still include the resonance contributions "0(4040), 0(4160) and 0(4415), which we 
subtract and include in the J/0-family resonance contribution. Finally from the T to 40 GeV we take weighted 
averages for the data sets JADE, TASSO, MD-1 and others as shown in Fig. 8 . Consistent results are obtained 
if we combine the data pointwise and add all systematic errors linearly. The results for aj) ad are presented in 
Table 3. The last two 
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Table 3a: Contributions to a dad • 10 10 

H' 


final state 

energy range (GeV) 

contribution (stat) (syst) 

rel. err. 

abs. err. 

P 

(0.28, 0.81) 

426.66 ( 5.61) (10.62) 

2.8% 

1.7% 

i0 

(0.42, 0.81) 

37.76 ( 0.45) ( 1.02) 

3.0% 

0.2% 

<t> 

(1.00, 1.04) 

38.55 ( 0.54) ( 0.89) 

2.7% 

0.1% 

JH> 


8.60 ( 0.41) ( 0.40) 

6.7% 

0.1% 

T 


0.10 ( 0.00) ( 0.01) 

6.7% 

0.0% 

hadrons 

(0.81, 1.40) 

112.85 ( 1.33) ( 5.49) 

5.0% 

0.8% 

hadrons 

(1.40, 3.10) 

56.43 ( 0.45) ( 7.22) 

12.8% 

1.0% 

hadrons 

(3.10, 3.60) 

4.47 ( 0.23) ( 0.86) 

19.9% 

0.1% 

hadrons 

(3.60, 9.46) 

14.06 ( 0.07) ( 0.90) 

6.5% 

0.1% 

hadrons 

(9.46, 40.0) 

2.70 ( 0.03) ( 0.13) 

4.9% 

0.0% 

perturbative 

(40.0, oo ) 

0.16 ( 0.00) ( 0.00) 

0.2% 

0.0% 

total 


702.35 ( 5.85) (14.09) 

2.2 % 

2.2% 


columns give the relative uncertainty (rel.err.) of the individual contribution and the absolute uncertainty 
(abs.err.) relative to the total result. The p-contribution always includes a contribution 2.08 (0.01)(0.05) 
calculated using chiral perturbation theory from the 27r-threshold to 318 MeV where data points start. For 
details we refer to Sec. [4.2j In spite of a major effort in the p region the result is almost the same as the one 
obtained using Kinoshita’s fit [Q which was used subsequently in |26[ |2?| ]. A previous analysis |^l[ 

gave the results 428.95(1.71)(12.81)[426.68(1.70)(12.74)] for the p contribution in the range (0.28, 0.81) GeV 
and 125.54(3.89)(11.61)[124.27(3.83)(11.46)] for the range (0.81,1.40) GeV. The values in square brackets are 
the ones obtained after correcting for the missing subtraction of the hadronic vacuum polarization. While the 
p contribution below 0.81 GeV remains unchanged the p-tail plus background up to 1.4 GeV turns out to be 
somewhat smaller. 


If we would treat the different experiments separately and then take weighted averages in common domains we 
would obtain 436.30(9.95)(8.02) for the p up to 0.81 GeV, a somewhat higher value. Details are illustrated by 
Tab. 3b. The high value is the result of a failure of 


Table 3b: p-contributions to a dad • 10 10 of different experiments 


Energy range: 

(0.32, 0.36) GeV 

(0.36, 0.40) GeV 

(0.40, 0.81) GeV 

DM1+ 

CMD 

OLYA 

6.213 (0.372)(0.214) 

9.225 (0.432)(0.197) 
11.533 (0.710)(0.231) 

409.739(14.450)( 9.331) 
414.176(11.051)( 8.284) 
472.399(24.364) (22.672) 

Averaged 

6.213 (0.372)(0.214) 

9.890 (0.373)(0.146) 

418.120( 9.939)( 7.661) 

Tab. 3a 

6.177 (0.370)(0.210) 

9.363 (0.582)(0.269) 

409.034( 5.572)(10.136) 


the trapezoidal rule when applied to the OLYA data. The reason is that the three lowest points have much 
larger energy separations than the other points at higher energies and a high weight due to the 1/s 2 behavior 
of the kernel. Therefore the first three points yield a large contribution to the integral. Actually, applying the 
trapezoidal rule here we overestimate the contribution because the low energy tail of the resonance is a strongly 
varying concave function. This problem is largely circumvented by combining data from different experiments 
before integration. Note that, while the integrals taken over the whole range agree quite well between different 
experiments, the individual contributions in the subdomains do not agree within errors. Apparently the integral 
is better defined than the local values. In other words, the p —> irn region is after all not as well established as 
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the experiments claim. There are many points where the weighted average yields S > 1 and the error must be 
enlarged. 

If we just take the collection of all points in a given energy region and add systematic errors linearly we obtain 
413.80(8.10)(23.81) for the range up to 0.81 GeV, with a much larger systematic error. Similarly, we find 
113.71(1.31)(12.31) for the contribution from 0.81 to 1.40 GeV. The final result would be 695.42(8.30)(27.52). 
We remind the reader that in the important region from 0.81 to 1.4 GeV data for the channels 7r + 7r _ 27r° and 
from ND lie substantially higher than from other experiments as can be seen in Fig. 4. A new 
experiment is required to resolve the current discrepancy. 

In Tab. 3c we give the results we obtain when using the resummed photon propagator in the diagram Fig. 1. 
This dressed contribution we denote by a dad *. We notice that the difference between the dressed and the 
undressed form is about 23, slightly larger than the uncertainty of 16 and of the same size as the interesting 
weak contribution a™ eak ~ 20 . For a T the results are collected in Tab. 4. The hadronic contribution to the 
electron anomaly is a kad * = (194.48 ± 1.69 ± 3.87) x 10 -14 . 


Table 3c: Contributions to a kad * • 10 10 

_ H' _ 


final state 

energy range (GeV) 

contribution (stat) (syst) 

P 

(0.28, 0.81) 

438.63 ( 5.76) (10.91) 

LU 

(0.42, 0.81) 

38.90 ( 0.47) ( 1.05) 


(1.00, 1.04) 

39.86 ( 0.56) ( 0.92) 



9.05 ( 0.43) ( 0.42) 

T 


0.11 ( 0.00) ( 0.01) 

hadrons 

(0.81, 1.40) 

116.84 ( 1.38) ( 5.69) 

hadrons 

(1.40, 3.10) 

58.92 ( 0.47) ( 7.55) 

hadrons 

(3.10, 3.60) 

4.71 ( 0.25) ( 0.90) 

hadrons 

(3.60, 9.46) 

14.91 ( 0.08) ( 0.96) 

hadrons 

(9.46, 40.0) 

2.94 ( 0.03) ( 0.14) 

perturbative 

(40.0, oo ) 

0.18 ( 0.00) ( 0.00) 

total 


725.04 ( 6.01) (14.57) 
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Table 4: Contributions to a!) ad * • 10 8 


final state 

energy range (GeV) 

contribution (stat) (syst) 

P 

(0.28, 0.81) 

141.35 (1.66) (3.42) 

LO 

(0.42, 0.81) 

15.06 (0.18) (0.41) 

<f> 

(1.00, 1.04) 

21.06 (0.29) (0.48) 

JH’ 


13.47 (0.64) (0.65) 

T 


0.26 (0.01) (0.01) 

hadrons 

(0.81, 1.40) 

58.23 (0.64) (3.11) 

hadrons 

(1.40, 3.10) 

58.84 (0.56) (7.94) 

hadrons 

(3.10, 3.60) 

7.02 (0.36) (1.35) 

hadrons 

(3.60, 9.46) 

28.04 (0.14) (1.80) 

hadrons 

(9.46, 40.0) 

7.36 (0.07) (0.35) 

perturbative 

(40.0, oo ) 

0.50 (0.00) (0.00) 

total 


351.18 (2.04) (9.51) 
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We now turn to the hadronic contribution to the shift in the fine structure constant. In Table 5a we list the 
contributions from different energy ranges. 


Table 5a: Contributions to Aa^ x 10 4 


final state 

energy range (GeV) 

contribution (stat) (syst) 

Ref. @ 

P 

(0.28, 0.81) 

26.08[26.23] (0.29) (0.62) 

26.07(0.10)(0.78) 

UJ 

(0.42, 0.81) 

2.93[ 2.96] (0.04) (0.08) 

3.43(0.35)(0.10) 

<t> 

(1.00, 1.04) 

5.08[ 5.15] (0.07) (0.12) 

5.27(0.24)(0.16) 

J/ 1 > 


11.34[11.93] (0.55) (0.61) 

10.16(1.34)(1.52) 

T 


1.18[ 1.27] (0.05) (0.06) 

1.17(0.04)(0.07) 

hadrons 

(0.81, 1.40) 

13.83[13.99] (0.15) (0.79) 

15.63(0.68)(1.73) 

hadrons 

(1.40, 3.10) 

27.62[28.23] (0.32) (4.01) 

27.95(0.60)(5.59) 

hadrons 

(3.10, 3.60) 

5.82[ 5.98] (0.30) (1.12) 

5.98(0.31)(1.15) 

hadrons 

(3.60, 9.46) 

50.60[50.50] (0.24) (3.33) 

50.27(0.51)(3.14) 

hadrons 

(9.46, 40.0) 

93.07 (0.86) (3.39) 

93.60(1.24)(2.91) 

perturbative 

(40.0, oo ) 

42.82 (0.00) (0.10) 

42.67(0.29)(0.59) 

total 


280.37[282.13](1.18) (6.43) 

282.21(2.19)(8.30) 


Table 5b: Shape dependence of resonance contributions 


final state 

energy range (GeV) 

a) 

b) 

c) 

d) 

LO 

(0.42, 0.81) 

3.138 

2.994 

2.952 

2.931 

<t> 

(1.00, 1.04) 

5.451 

5.068 

5.069 

5.083 

J/i> 


11.380 

11.342 

11.338 

11.338 

T 


1.182 

1.178 

1.178 

1.178 


a) narrow width approximation 

b) non-relativistic 

c) relativistic constant width 

d) relativistic s-dependent width 


The resonance contributions were evaluated as discussed in Sec. 4.1. Note that Breit-Wigner resonances may 
be treated non-relativistically, relativistically and with different off-resonance behavior. Results obtained for 
different resonance shapes are listed in Tab. 5b. We note that the narrow width approximation gives generally 
larger values than the Breit-Wigner parametrizations in either the non-relativistic form, the relativistic form 
with constant width or the relativistic form with s-dependent width. As expected, the deviations obtained from 
using different types of Breit-Wigner parametrizations are within the experimental uncertainties. 


The region between the J/ip and the T has been split into two subdomains. From the J/x[> to 3.6 GeV only 
Mark I data are available. Above 3.6 GeV up to the T one may think about skipping the Mark I data, as 
mentioned before. Because of the resonances included in the data we integrate the PLUTO, DASP and Mark I 
data separately in the region from 3.6 to 5.2 GeV and combine the results after integration. 


The analysis shows that the results are affected in a minor way if we include the Mark I data. Taking the 
weighted average of the two values obtained from the Mark I data on the one hand and the other data on the 
other hand we find 50.79(0.20)(3.20). If we combine results from individual experiments and take weighted 
averages for overlapping domains we get 50.69[51.20](0.30)(2.74) if we exclude [include] the Mark I data. These 
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checks show that the trapezoidal rule works consistently and there is no serious problem of properly weighting 
the data from different experiments according to their uncertainties^. The values obtained compare with the 
ones from a previous analysis |^7j and are all found to be consistent. 

We have checked the non-resonant part above the J/ij), which contributes a major part to A a. against the 3- 
loop perturbative QCD prediction, using the LEP value a s (M|) = 0.126 ±0.005 (23) for the strong interaction 
constant as an input. We obtain 


Range 


Data 

Ref. H 

Eq. (|) 

5.00 - 9.46 

GeV 

32.63 

35.78 

34.97 

12.0-40.0 

GeV 

79.22 

78.23 

77.64 


The two QCD results differ by a different treatment of mass effects. In Eq. (||) just the lowest order threshold 
factor is used, while Ref. |Q also takes into account charm and bottom mass effects in the higher order terms 
in an expansion in m 2 / s. 

We have also varied the high energy cut energy from E cut =40 to 30, 20 and 12 GeV and found stable results: 
0.0280 ± 0.0007, 0.0280 ± 0.0006, 0.0280 ± 0.0006 and 0.0279 ± 0.0006. We note that data and the perturbative 
prediction in average fit fairly well above the resonance regions, as may be seen in Figs. 7 and 8 . In addition we 
observe that the onset of the 7 — Z mixing is well under control after the subtraction of the 7 — Z interference 
term which we described at the end of Sec. |3] This subtraction was applied already in previous work |^ 6 |, |27| |. 
Doubtless, the LEP experiments have dramatically improved our confidence in perturbative QCD and we may 
well use a much lower E cut , such as for example 12 GeV, and use perturbative QCD also in the range from e.g. 
5 GeV, up to the T threshold. As a result we find 0.0282 ± 0.0005 and hence a slightly larger result with a 
smaller error. However, if we move the cut slightly from 5 GeV to 4.5 GeV we obtain 0.0280 ± 0.0005 which 
shows that the better accuracy is delusive as the central value depends substantially on the cuff]. A glance at 
Fig. 7a shows that the PLUTO data points in the region of the -0(4040), 0(4160) and -0(4415) resonances are 
lower than “perturbative QCD plus the Breit-Wigner resonances”. Note that the resonances depicted in Fig. 7a 
have been scaled down by ( a/a(s )) 2 in order to have the proper normalization for R. If one applies the missing 
(1 ± 2A ai){a/a{s)) 2 correction to the Mark I data they agree much better than the PLUTO data with the 
above “prediction”. 

As a main result we obtain 


Aa£2j = 0.0280[0.0282] ± 0.0007 (21) 

for Mz = 91.1888 GeV consistent with the 1991 update [^?J. In brackets we also give the value without 
any ( a/a(s )) 2 rescaling of data. These results confirm the assumption made in previous work, namely, that 
corrections which account for the missing subtractions of the vacuum polarization contributions are small and 
well within errors. In fact usually the “missing corrections” were estimated and included as a part of the 
systematic error. This small correction was not taken into account in j27j. 

3 This is in contrast to claims in Ref. j|^] that the trapezoidal rule does not allow us to combine data from 
different experiments in a reasonable way and automatically has the effect of weighting all inputs equally. Note 
that in previous analysis 0 0 -H properly weighted averages of data from different experiments were used 
unless data from different experiments had very similar errors in which case they were treated like points from 
one experiment. As systematic errors for a given “experiment” are added linearly the simplified treatment of 
equal weighting of data has primarily the effect of yielding a too conservative estimate for the systematic error. 
In Ref. (28) data are fitted to smooth functions before integration and systematically lower results are obtained 
(see also the discussion in Ref (bl]]). The author of Ref. (^ 8 | believes more in integration of his fits than in 
trapezoidal integration, without giving any actual arguments; moreover he does not mention the fact that his 
fitting method has to rely on some rather arbitrary assumptions about the sources of systematic effects and/or 
correlations. 

4 Our findings do not confirm those of Ref. ^9|. These authors essentially use perturbative QCD with the 
world average strong coupling a s = 0.118 ±0.007 for E > 3 GeV, plus the resonances, instead of the data. Data, 
in so far as they are used, are rescaled to meet the result of perturbation theory. Their results lie systematically 
lower than ours. The third order QCD prediction is assumed to be exact down to 3 GeV. 
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A lower central value than in |24|, [26] is mainly due to the new more precise data from CB [IOC] and MD-1 [101]. 
The Mark I data dominated in most of the earlier estimates. The replacement of the Mark I data by the CB 
data in the common energy range has been presented in the update |27j some time ago. Undressing of the 
resonance contributions and some supplementary subtractions of hadronic vacuum polarization contributions 
lead to a correction of —0.000176. We now employ the value for a s from LEP. In Refs. |24|, ^6|, ^7j the a s value 
used was the one obtained by the PEP/PETRA experiments. 

The uncertainty obtained in this analysis is smaller mainly for the following reasons. First, we are using more 
complete data in the range from 0.81 to 3.1 GeV. In particular the DM2 data J74| in the range from 1.35 to 
2.3 GeV, which have smaller uncertainties, were not used in previous estimates. Second, previously in [f26| , |27[ 
an overall 20% systematic error was assumed in this range, which corresponded to a typical systematic error 
reported by the individual experiments. Furthermore, in these references, a conservative 3.5% error, which 
was the accuracy of the PEP/PETRA data, was assumed for the perturbative tail. This more conservative 
treatment of errors in Ref. [£6| lead to an increase to 0.0009 from 0.0007 which was previously obtained in 
Ref. [Q. The more complete collection of data allows us to better justify the more precise result now. Although 
the replacement of the Mark I data by the CB data helps in reducing the uncertainty it is not the main reason 
as can be checked in Tab. 5a. The higher accuracy of the resonance parameters also helped in reducing the 
uncertainty. 

We finally summarize in Tab. 5c the uncertainties obtained for the different contributions. 


Table 5c: “Distribution” of uncertainties 



A«iad >< IQ" 

rel. err. 

abs.err. 

Resonances: 

46.61 (1.08) 

2.3 % 

0.4 % 

Background: 




E < Mj 

41.45 (4.11) 

9.9 % 

1.5 % 

Mj^ < E < 3.6 GeV 

5.82 (1.16) 

19.9 % 

0.4 % 

3.6 GeV <E <M r 

50.60 (3.34) 

6.6 % 

1.2 % 

M r < E < 40 GeV 

93.07 (3.50) 

3.8 % 

1.2 % 

E < 40 GeV data 

237.55 (6.54) 

2.8 % 

2.3 % 

40 GeV < E QCD 

42.82 (0.10) 

0.2 % 

0.0 % 

total 

280.37 (6.54) 

2.3 % 

(*) 

(4.92) 

(1.8 % ) 


The last line (*) of Tab. 5c gives the uncertainty one would get if the experimental errors on R(s) would 
be reduced to 5% in all the regions which exhibit uncertainties larger than that. The last column gives the 
uncertainty relative to the total result. The table clearly tells us that the background contributions require a 
new scan for all energies up to about 12 GeV. For higher energies one may rely more on perturbative QCD. 
It seems unlikely that a substantial improvement will be possible in the foreseeable future. Physics will need a 
global update of many experiments in order to be prepared for the next level of precision physics. 


7. Summary and outlook 


The question at the beginning of this investigation was whether one could improve the previous estimates 
in the present situation and how things could develop in future. While there has been very little truly new 
experimental results, some groups have published updated results which were available before as preliminary 
data only. Examples are the ND results {54|, the DM2 results |6^, |74|, [75], [fO and the Crystal Ball results 10C], 
where the latter have already been used in [E7|. New results from VEPP-4 [101 have been included as well. In 


addition we have made an effort to collect as much as possible all the available data. It should be noted that 
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for example the Durham-RAL High Energy Physics database utilized in Ref. (29J is incomplete and in many 
cases contains preliminary data only while the final results are missing. In our analysis we also used data which 


have not yet been included in the collection [ 1111 of e + e^ data which was published recently. Another issue was 
to check what corrections were applied to the published data, which required a careful reading of the original 
papers. 

Additional motivations for performing this update were the following: the recent issues of the Review of Particle 
Properties ]4l|| had some improvements on the resonance parameters and there was progress in calculating R(s) 
at high energies in perturbative QCD 


p3L 231. For the muon anomaly the study of the low energy end by 


means of chiral perturbation theory |48| allowed us to reduce potential model-dependences of earlier approaches. 

That all these efforts mentioned lead to minor changes of the results was to be expected. Nevertheless we think 
such an update was necessary and it will not be soon that a more precise estimate will be possible. 

Refinements and improvements were proposed recently in R efs. [^8l and [pij. More theoretical attempts to 
calculate the photon vacuum polarization may be found in [ 13_2 , 113 , |l!4j , 115 , 116 1 to mention only a few. 
We all know about the difficulties to make accurate reliable predictions in strong interaction physics. Here we 
have tried to estimate hadronic vacuum polarization effects in a model-independent way exploiting as much as 
possible the existing data, and we find results consistent with known estimates. 


For the muon anomaly we propose to use the RG improved value which is ler higher than the bare one usually 

and assuming that the subleading hadronic contributions and their 


considered. With our value for a!l ad * 


uncertainties are as given by Kinoshita et al. fllq] , we obtain 


a^ he = (11659210 ± 16) x 10 


-10 


and therefore at present we have 

a® xp - aj) he = (20 ± 86) x 10" 10 . 

The forthcoming Brookhaven experiment is expected to reduce this uncertainty to ±17 x 10 -1 ° . Recently there 
is an ongoing discussion about the true size of the hadronic light by-light scattering contribution |3l[] and 
therefore the value and the uncertainty of the theoretical prediction still may change. 

Finally, by adding the hadronic contribution, the leptonic contribution Acq = 0.031421 and a top contribution 
Acq op = —0.000061 (obtained by integrating Eq. (|[) with running parameters numerically) we find the total 
shift by the fermionsjj 

Aa = 0.05940 ±0.00065 

and hence for the effective fine structure constant 


a(Mf)- 1 = 128.896 ±0.090 . 


What about the future? Progress in decreasing the hadronic uncertainties can be expected from the experiments 
at the general purpose CMD-2 detector at the Novosibirsk VEPP-2M collider, which has an ambitious goal of 
measuring <Jtot(e + e~ —> hadrons) with an accuracy better than 1%. The experiment is in progress and during 
the 1995 runs it is planned to study the energy range from the threshold of two pion production up to the <fi 
meson with a 10 MeV step and achieve a statistical accuracy of 3% in each point |ll8|] . Further work will be 
needed to understand the detector performance as well as low energy radiative corrections at such a level that 
systematic uncertainties will be understood with the desired accuracy. After that the energy range from the <fi 
meson up to the maximum attainable energy of 1.4 GeV will be studied. Here, as discussed above, additional 
difficulties can arise for the high precision cross section measurements because of the intermediate mechanism 
uncertainty. Let us assume that a systematic uncertainty of 1% will be achieved for the p,u and <p mesons and 
that of 2% for the hadron continuum between 0.81 and 1.40 GeV. The statistical error in the integrals can be 
neglected and a resulting systematic uncertainty in the muon anomaly decreases from 15 x 10~ 10 to 9 x 10 -1 ° 
of which 8 x 10 -1 ° come from the contribution of the hadron continuum between 1.4 and 3.1 GeV. 


Further progress can be expected from the future generation experiments at DA4>NE at Frascati and the 
upgraded VEPP-2M. Plans at DA4>NE are to perform a scan with 100 values from 0.28 to 1.5 GeV c.m. energy 


5 We do not include the W boson contribution here because it is gauge dependent, and splitting off a gauge 


invariant part and combining the reminder with the photon vertex is not unique 1117 
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with a precision which will allow to reduce the uncertainty from the hadronic vacuum polarization contribution 
to a M to about 0.3% ]119[ . This will reduce the uncertainty for the region below 1.4 GeV to 1.5 x 10 -10 and the 
remaining total hadronic vacuum polarization uncertainty in the prediction of a^ would be 8 x 10~ 10 , which 
again is completely dominated by the contribution from 1.4 to 3.1 GeV. The DA4>NE measurement will be 
possible earliest by end of 1997 but then can be done within a few days. Obviously, the VEPP-2M and the 
DA$NE measurements will be crucial for the physics we will be able to learn from the planned Brookhaven 
Experiment 821 [ p~2| . 

No improvement is in sight at higher energies. The energy region from 1.4 to 3.1 GeV and in addition that 
of higher energies from 3.1 up to about 12 to 40 GeV, depending on how much one is accepting to rely on 
perturbative QCD, will give a dominant contribution to the uncertainty of the fine structure constant shift. 

One can conclude that a real breakthrough in improving the precision of the hadronic vacuum polarization will 
require dedicated efforts in high precision R-measurements in a wide energy range. 
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Appendix: “Dressed” versus “undressed” quantities 

Equations (jli) and (0) can be derived from the convergent dispersion relation 

Imn^(s') 


Re An' (s) = Re lT(s) - lT(0) = - Re 


ds' 


s'(s' — s — is) 


where n; is the transversal part of the current-current matrix element of the conserved electromagnetic current 
iW q) =i j d A x e il * x < 0| T* j 7 (z) jj(0) |0 >= (<?V - q»q v )il\(q 2 ) . 

This self-energy function exists in the limit where the electromagnetic interaction is switched off and for the 
hadronic current, in principle, it is determined by QCD. The irreducible photon self-energy is given by Ii;(s) = 
e 2 n;(s). It is well known that the running of a = e 2 /47r (see Sec. 3) may be understood as a consequence of 
the Dyson summation of the irreducible photon self-energy 


1 


d (1 + R ea n; W )^ R , T -^- J 


e 2 (») 


which is equivalent to a RG improvement of the perturbation expansion. For the full photon propagator we 
have 


Re ■ 


1 


1 


1 


1 


e 2 (s) 


Irn II' (s) 


1-An;(s) i-ReAn;(s) 1 + f _ 

7 7 1 + 11—Re AIT^(s) 

for the real part and 

1 _ Imli; (s) 1 


Y 1 — Re An;(s) 


Im 


— im n;(s) 


e\s) 


l-AIE(s) (1 — Re Ant(s)) 2 i , « 

7 7 1 + 11 —Re AIli,(s) 1 

for the imaginary part. 

These relations show that in the dispersion relation for the irreducible self-energy one has to use “undressed” 
(with respect to vacuum polarization effects) quantities and not the physical ones. Thus in the optical theorem 
(unitarity) one has to use the undressed total cross section 

- Im n; (s) = • 
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Up to subleading higher order effects the resummed photon propagator may be obtained directly from a dis¬ 
persion relation by using the dressed physical quantities under the dispersion integral. 
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Fig. 3: The pion form factor |F,r(s)| 2 for 0.32 to 0.81 GeV. In this and the following 

figures only the statistical error bars are shown. 
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Fig. 4: Averaged value of R for 0.81 to 1.4 GeV calculated from the exclusive channels 

using all available data. For comparison the R values compiled in Ref. |Q 
are shown. 


Fig. 5: R{s ) in the range 1.4 to 2.3 GeV calculated from the available R{n > 2) 

compilations with the two-body channels added up. 
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Fig. 6: 


The R data from the various experiments in the energy range 3.1 to 9.6 GeV. 


Fig. 8: The R data from the various experiments in the energy range 9.6 to 40 GeV. 

The two perturbation theory results are the 0{a 3 s ) predictions Eq. (||) and 
Ref. jlH with the LEP value a s (M§) = 0.126 ± 0.005 as input. 
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Fig. 7: The R data from the various experiments in the energy range 3.1 to 9.6 GeV. 

The “weighted mean” does not include the Mark I data above 3.6 GeV. The 
two perturbation theory results are the O(af) predictions Eq. ([!]) and Ref. |Q 
with the LEP value a s (M= 0.126 ± 0.005 as input. 
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